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ABSTRACT: The Bragg regime shifts when conventional materials in a multilayer 
distributed Bragg reflector (DBR) are replaced by artificial materials with the so- 
called negative index of refraction. This provides an avenue for characterizing the 
latter class of materials. 
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Artificial materials with the so-called negative index of refraction at ~ 10 GHz fre- 
quency have garnered much recent interest These isotropic materials, with 
supposedly negative real permittivity and negative real permeability, were the subject 
of a mid-1960s paper of Veselago [^, who predicted certain unusual electromagnetic 
properties and coined the unclear term left-handed materials for them. Handed these 
materials are not P] , issue can be taken on the isotropy of their first samples , and 
strictly non-dissipative no materials can be P]. Yet the only available experimental 
result § conclusively shows that these materials are different from their conventional 
counterparts (i.e., those with positive real permittivity and positive real permeabil- 
ity). As better realizations are probably in the wings, in colloquial terms, the business 
of macroscopic electromagnetics can no longer be as usual. 

The manufacturing process of the subject artificial materials delivers samples in the 
form of slabs. Even a cursory perusal of key papers ll|]-[0] reveals that characteristic 
observable properties of these materials supposedly are manifested most clearly when 
plane waves are obliquely incident on planar interfaces with their conventional coun- 
terparts. Indeed, the only experimental confirmation of their unusual characteristics 
exploited oblique incidence on a planar interface with air [0. In contrast, we present 
here a configuration wherein the incidence is normal and yet the subject artificial 
materials can be easily distinguished from their conventional counterparts as follows. 
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Distributed Bragg reflectors (DBRs) are commonly used planar devices in optics |]T0|- 



IJ]. A multilayer DBR is a stack of layers of two alternating materials with different 
indexes of refraction as well as low absorption, exhibiting very high reflectance in 
the so-called Bragg regime. If P denotes the thickness of a unit cell comprising two 
adjacent layers, and n is some effective index of refraction, then the Bragg regime is 
located about the (free-space) wavelength [|ll|]-[|l^ 

Af = 2Pn. (1) 

(Bragg regimes are also possible around integer submultiples of this Af*", but it must 
be borne in mind that indexes of refraction are wavelength-dependent, in general.) 
The parameter n depends on the indexes of refraction and the volumetric proportion 
of the two constituent materials. 

Suppose the region < 2; < L is occupied by a multilayer DBR, as shown in Figure 
1. The multilayer DBR comprises N unit cells, each made of two layers labeled a and 
b, with relative permittivities ea^t, and relative permeabilities fia,b- With P = L/N 
as the thickness of the unit cell, the thickness of layer a is equal to qP, < q < 1. 
A plane wave is normally incident on the DBR from the vacuous half-space z < 0, 
with Ao denoting its wavelength. Therefore, a reflected plane wave also exists in the 
half-space z < 0, and a transmitted plane wave in the vacuous half-space z > L. 
The corresponding electric field phasors are given by 

(aexriiikoz) + rexp(— zA;„z) , z < 

- , (2) 
texp [ik(-,{z — L)] , z > L 

where kg = 27r/Ao is the free-space wavenumber; a, r and t are the amplitudes of 
the incident, reflected and transmitted plane waves, respectively; while (u^., Uy, u^) is 



Gerardin and Lakhtakia: Negative index of refraction and DBRs 



4 



the triad of cartesian unit vectors. An exp{—iujt) time-dependence is implicit, where 
uj = /co/(eo/io)^/^ is the angular frequency, while eo and /Iq are the permittivity and 
the permeability of free space, respectively. 

The amplitudes r and t must be determined in terms of a. This is best done by 



setting up the 2x2 matrix equation 15 
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where rf^ = (/io/co)^/^ is the intrinsic impedance of free space. This equation has to 
be numerically solved, which we did. The principle of conservation of energy requires 
that + < with the equality holding only if the both constitutent materials 
in the DBR are non-dissipative at the particular value of Aq. Our algorithm satisfied 
the conservation principle. 



3 NUMERICAL RESULTS 

Figure 2 shows the computed reflectance |r/ap as a function of Xq/P for three values 
of q when e„ = 4(1 + zO.OOl), fia = 1-02(1 + zO.OOl), = 1, /i^ = 1 and = 20. 
The Bragg regime is clearly observable via the rectangular feature with an almost 
fiat top and two vertical sides in all three plots. The full-width-at-half-maximum 
bandwidth AAq of the Bragg regime is typically < 1.25 P for all q. Predictably, the 
Bragg feature vanishes for g = and q = I- As q increases, so does \f^/P (at the 
center of the Bragg feature in each plot); which amounts to an increase in n, as shown 
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in Table 1. 

Ideal Bragg features do not emerge for all values oi q E [0.2, 0.5], when = 
4(-l + i0.001), Ha = 1.02(-l + i0.001), = 1, /Xb = 1 and = 20. Thus, the Bragg 
feature for q = 0.5 is not flat-topped in Figure 3, although it is well-developed for 
q — 0.6 and q — 0.7. Calculated values of n as functions of q are shown in Table 1. 

We note from the presented and related results that n > 1 for all values of g e (0, 1) 
when Re[ea, > 0. In contrast, < n < 1 for ^ ~ 0.67 when Re[ea, < 0. The 
reduction of n below the unit index of refraction of material b could suggest that the 
real part of the index of refraction is negative for the subject artificial materials, but 
that suggestion does not appear to be supported by the values of n > 1 for g ~ 0.67 
when Re[ea,//(j] < 0. Anyhow, in conjunction with Figure 3, Table 1 confirms that 
wave-material interaction in the subject artificial materials is intrinsically different 
from that in their conventional counterparts. 

Our results also show that the Bragg regime would shift to shorter wavelengths, 
if a conventional dielectric/magnetic constituent of a multilayer DBR were to be 
replaced by its analog of the subject variety. Consequently, measurements of n would 
illuminate the issue of the negative index of refraction, and could also help in the 
characterization of the subject artificial materials. At the same time, multilayer DBRs 
made with the subject artificial materials could be useful in wavelength regimes that 
are inaccessible with DBRs made with only conventional materials. 
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Table 1 Values oi n — Xf^/2P computed for different values oi q; N — 20, €b — 1 
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e„ = 4(l + i0.001) 


e„ = 4(-l + i0.001) 




fXa = 1.02(1+^0.001) 


fla = 1.02(-l+i0.001) 


0.1 


1.14 


0.69 


0.2 


1.26 


0.37t 


0.3 


1.38 


0.09t 


0.4 


1.47 


0.23t 


0.5 


1.58 


0.56^ 


0.6 


1.68 


0.82 


0.7 


1.78 


1.09 


0.8 


1.87 


1.37 


0.9 


1.96t 


1.68^ 



^ Bragg feature has a curved top — exemplified in the top plot of Figure 3 — which 
begins to flatten as N increases. There exists also a certain degree of arbitrariness in 
the identification of the Bragg feature for some values of q between ~ 0.25 and ~0.42 
when Re[ea] = —4 and Re[/ia] = —1.02. 
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Figure 1: Schematic of the boundary value problem. The distributed Bragg reflector 
comprises unit cells, each of thickness P and containing one layer each of two 
different materials labeled a and b. 
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Figure 2: Reflectance |r/ap of a DBR as a function of Aq/P for different values of g; 
AT = 20, ea = 4(1 + iO.OOl), /x^ = 1.02(1 + zO.OOl), = 1 and /Xb = 1. 
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